Fast magnetization switching of Stoner particles: A nonlinear dynamics picture 
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The magnetization reversal of Stoner particles is investigated from the point of view of nonlinear 
dynamics within the Landau-Lifshitz-Gilbert formulation. The following results are obtained. 1) 
We clarify that the so-called Stoner- Wohlfarth (SW) limit becomes exact when damping constant 
is infinitely large. Under the limit, the magnetization moves along the steepest energy descent 
path. The minimal switching field is the one at which there is only one stable fixed point in the 
system. 2) For a given magnetic anisotropy, there is a critical value for the damping constant, above 
which the minimal switching field is the same as that of the SW-limit. 3) We illustrate how fixed 
points and their basins change under a field along different directions. This change explains well 
why a non-parallel field gives a smaller minimal switching field and a short switching time. 4) The 
field of a ballistic magnetization reversal should be along certain direction window in the presence 
of energy dissipation. The width of the window depends on both of the damping constant and 
the magnetic anisotropy. The upper and lower bounds of the direction window increase with the 
damping constant. The window width oscillates with the damping constant for a given magnetic 
anisotropy. It is zero for both zero and infinite damping. Thus, the perpendicular field configuration 
widely employed in the current experiments is not the best one since the damping constant in a real 
system is far from zero. 

PACS numbers: 75.60.Jk, 75.75.+a, 05.45.-a 



I. INTRODUCTION 

Magnetic data storage is one of the important com- 
ponents of modern computers. Data input and output 
involve switching the magnetization of magnetic stor- 
age cells (magnetization reversal). The typical switch- 
ing time with currently used technology is of order of 
nanosecond. If one wants to have a faster computer 
(modern electronic computers are working at a clock 
speed of order of GHz), the conventional magnetization 
reversal method shall soon (the clock speed is doubled 
every year in the past) become a bottleneck. Thus fast 
magnetization switching shall be of great importance for 
the future development of high speed information indus- 
try. 

The reversal of a magnetization can be achieved in 
many different ways, and it is a very complicated issue 0. 
For example, in a bulk material, the magnetization rever- 
sal can go through bucking and curling modes, or nucle- 
ation and domain formation. The recent advance of tech- 
nology allows us to fabricate the magnetic nano-particles 
that are believed to be useful for high density informa- 
tion storage0 El For a magnetic nano-particle, 
the magnetic moments of all atoms are normally aligned 
in the same direction, creating a so-called single mag- 
netic domain. Such a nano-particle is usually called a 
Stoner- Wohlfarth or Stoner particle. The understanding 
of magnetization reversal of a single magnetic domain 
should be relatively simple in comparison with that in a 
bulk system, but important in nano-technologies0. 

There are two challenging issues about magnetization 
reversal. One is how to have a short reversal time, and 
the other is how to make the switching field to be small. 



The conventional magnetization reversal technique is to 
apply a magnetic field antiparallel to the initial magne- 
tization. A large enough field can drive the initial state 
out of local minimum and at the same time make the 
target state to be the global minimum. Thus the sys- 
tem can roll down to the target state through ringing 
effect 0, 0, 0- For the issue of minimal switching field, 
the classical result, called Stoner- Wohlfarth (SW) limit, 
was given by Stoner and and Wohlfarth [Io|. 

Recently picoseconds ma gne tization switching has 
been observed in experiments |lllll2l | by using pulse mag- 
netic fields. Unlike the conventional method, the mag- 
netic field is applied in a perpendicular direction such 
that the magnetization undergoes a precession. This ap- 
proach has also received many theoretical attentions [UJ, 
llil Hil H(| . Numerical investigations [l^ showed that 
the switching time can be substantially reduced because 
ringing effect can be suppressed so that the magnetiza- 
tion will move along a so-called ballistic trajectory [li). 
The precessional magnetization reversal provides not 
only a shorter time but also lower switching field (well 
below the SW-limit), as found in the early numerical 
calculations 17]. In the absence of energy dissipation, 
precessional magnetization switching can also be inves- 
tigated analytically. Analytical results for the minimal 
field were obtained by Porter Recently Xiao and 
Niu have also studied the minimal field required in the 
precessional magnetization switching in a conservative 
system 19]. Their results were based on fast switching 
such that the energy dissipation can be neglected and the 
magnetic system can be regarded as a conservative sys- 
tem. They used the phase plane to reveal the properties 
of magnetization reversal. It was shown that precessional 
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magnetic switching occurs when localized trajectories in 
phase plane become delocalized. 

Although it is reasonable to approximate a Stoner par- 
ticle as a conservative system in a short time if the damp- 
ing (dissipation) during its dynamical motion is small. A 
real magnetic particle is not a truly conservative system 
when the damping is large or if one is interested in the 
long time behave. Dissipation should be taken into ac- 
count and a Stoner particle should be treated as a non- 
conservative system. In this paper we re-examine the 
magnetization reversal as a nonlinear dynamical system 
in the presence of dissipation. The magnetic dynamics 
of a single domain magnetic particle can be described by 
the evolution trajectories in the phase plane. One can 
use the general concepts of nonlinear dynamics not only 
to understand all results from previous studies, but also 
to see the validity conditions of some of these results such 
as the SW-limit. The paper is organized as follows. In 
Sec. II we first reformulate the magnetization reversal 
in terms of nonlinear dynamics concepts, such as attrac- 
tors and phase flow. Previous results are re-interpreted 
in such language. The Landau-Lifshitz-Gilbert (LLG) 
equation is introduced. Our main results are presented 
in Sec. III. The conditions under which the SW limit 
is valid are given. For a given magnetic anisotropy, we 
shall show that there is a critical damping constant above 
which the minimal switching field is the same as that of 
SW-limit. The reason and meaning of this critical damp- 
ing constant are also given. The change of fixed points 
and their basins under a magnetic field is investigated. 
We show that the field corresponding to the ballistic re- 
versal changes from the perpendicular direction for a con- 
servative system to a direction window in the presence of 
energy dissipation. The conclusion is given in Section IV. 

II. A NONLINEAR DYNAMICS PICTURE OF 
MAGNETIZATION REVERSAL 

II. 1 Attractors, phase flow and magnetization 
reversal 

Previous results on the magnetization reversal of 
Stoner particles can be conveniently described in the ter- 
minology of nonlinear dynamics. The phase space related 
to the magnetization is a two dimensional (2D) plane be- 
cause all atomic magnetic-moments are aligned in the 
same direction, and the magnetization can rotate under 
an external and/or an internal effective magnetic field 
that can exert a torque on the spin. The polar angle 9 
and the azimuthal angle (f>, shown in (Fig.^,), can fully 
determine a magnetization to. In the 9 — <j) plane, each 
point corresponds to a particular state of the magnetiza- 
tion. A state will in general evolute to new states due to 
its dynamics. Its motion can be described by a trajectory 
in the phase plane, called phase flow. The phase flow for 



a dissipative system ends to a few types of destiny (at- 
tractors), including fixed points, limit cycles, or strange 
attractors. They correspond to stable states, periodic, 
aperiodic and chaotic motions [20|. In a 2D phase plane, 
however, strange attractor solution is not allowed. 

The only attractor related to the magnetization rever- 
sal of Stoner particles is fixed points. The magnetization 
reversal problem is as follows: Before applying an exter- 
nal magnetic field, there are two stable fixed points (de- 
noted by A and B in Fig.^ 3 ), corresponding to magneti- 
zations, say Too (point A) and —too (point B), along its 
easy axis. The phase plane can be divided into two parts, 
called basins of attractors. One is around A, and the 
other around B, denoted by shadowed areas in Fig. ^d. 
The system in basin A(B) will end up at state A(B). Ini- 
tially, the magnetization is too, and the goal is to apply a 
small external field to switch the magnetization to —too 
fast. 




FIG. 1: (a) The magnetization rh can be uniquely deter- 
mined by angles 9 and <j>. z-axis is assumed to be along the 
total magnetic field h t . —rh x h t determines the precession 
direction, and — rh x (rh x ht) decides the dissipation direc- 
tion, (b) The 9 — cj> phase plane for the magnetization of a 
Stoner particle. Point A and point B represent the initial and 
the target state, respectively. Two shadowed areas denote 
schematically basins of two stable fixed points A and B. The 
solid curve LI and dashed curve L2 illustrate two different 
phase flow connected A and B. 



II. 2 Damping and non-damping magnetization 
reversals 

The conventional magnetization switching is based on 
a damping mechanism through, typically, the spin- lattice 
relaxation. From the point of view of nonlinear dynam- 
ics, the idea behind the method is to construct the ex- 
ternal magnetic field in such a way that the target state 
is the only stable fixed point. In another word, basin A 
(Fig. nj)) is reduced to zero and the whole (except proba- 
bly a few isolated points) phase plane is the target state 
basin (Basin B). The minimal reversal field (SW-limit) 
is the one at which basin A shrinks to a point. Since the 
initial and the target states have very large energy differ- 
ence, the extra energy must be dumped into the lattice 
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during a spiral motion before the system reaches the final 
state. The system first spirals out of A, and then spirals 
toward B, denoted by phase flow L2 in Fig.^i. This spi- 
ral motion is often referred [l^ to as ringing effect. The 
reversal time is largely determined by the effectiveness of 
energy dissipation. 

On the other hand, the target state does not need to 
be the only fixed point in the recent fast magnetization 
switching. In fact, it does not even need to be a fixed 
point. In the precessional magnetization reversal, one 
applies a short magnetic field pulse such that both initial 
and final states are not fixed points, and system will start 
to flow in the phase plane. In order to switch the magne- 
tization, one needs to let the system to reach the basin of 
the final state (Basin B) such that the system will flow 
to the target state after the pulse field is switched off. 
Ideally one wants both initial and target states on its 
precessional path. This is a non-damping method, and 
the reversal time does not rely on the spin relaxation 
time. There are several ways with different control pre- 
cisions to move the system to the desired state. One way 
is to apply a perpendicular pulse field to 'kick' the sys- 
tem to basin B. In comparison with conventional method, 
the spiral motion out of the initial state is replaced by a 
ballistic jl6| motion. However, the system relies on ring- 
ing effect to reach the final state. It was shown [l3^ that 
the switching time can be reduced substantially, but it 
is still hundreds of picoseconds for a normal magnetic 
particle due to the ringing effect in the last stage of mag- 
netization reversal. 

The proposal of Xiao and Niu^| is a fully ballistic 
magnetization reversal scheme. Neglecting the energy 
dissipation during the motion of a Stoner particle in a 
magnetic field, the magnetization must move on an equal 
potential line. Thus the idea is to apply an external mag- 
netic field with a proper strength and in a right direction 
such that both the initial and the target states have the 
same energy and are on the same phase flow trajectory as 
schematically illustrated by the solid line LI connected 
points A and B with an arrow in Fig.^. Without damp- 
ing, there is no extra energy to dissipate in this new ap- 
proach, thus the system move from the initial state to 
the final one in a ballistic way instead of in a ringing 
mode. The typical time for a precession of 180° in a field 
of teslas is order of picoseconds for usual magnetic mate- 
rials so a picoseconds magnetic field pulse is required in 
this method. As soon as the system arrives at the target 
state, one needs to switch off the external field, so the 
target state becomes a fixed point of the system again. 
This approach thus requires a precise control of the pulse 
duration. 

From computational point of view, the magnetization 
reversal time can be evaluated as soon as the phase flow 
connected the initial and the target states is found. The 
reversal time is given by the length dl of phase flow line 



divided by the phase velocity y 9 2 + ip 2 which is deter- 
mined by the system dynamics, 



t 



dl( 



(1) 



In the language of nonlinear dynamics, an external field 
modifies the dynamics by changing the phase velocity 
field. This velocity field is in general a continuous func- 
tion of the external field. A phase flow between the ini- 
tial and target states could only been set up when the 
external field is strong enough because the initial and fi- 
nal states are two stable fixed points with equally large 
basins at the beginning. The minimal switching field is 
the critical one at which such a flow is created. 



II. 3 The Landau-Lifshitz-Gilbert equation 

The magnetization dynamics of a Stoner parti- 
cle is governed by the Landau-Lifshitz-Gilbert (LLG) 
eauation|2l). 



which can also be written as 

(i + a 2 )^r = -\i\My.Ht- 



a - dM 

W Mx W 



M. 



M x (M x H t ). (2) 



Here I7I = 2.21 x 10 5 (rad/s)/(A/m) is the gyromagnetic 
ratio, M s is the saturated magnetization of the particle, 
and a is the phenomenological dimensionless damping 
constant. The typical experimental values of q|ll| ranges 
from 0.037 to 0.22 for different Co films. The equation 
describes the motion of magnetization vector M under an 
applied magnetic field H in the presence of an arbitrary 
magnetic anisotropy energy density function W{M,H). 
The total field H t = -V a W (M , H) / fi = Heff + H, 
where H e f / denotes the internal effective field due to the 
magnetic anisotropy. /io = 4ir x 10~ 7 N/A 2 is the vacuum 
magnetic permeability. 

It is convenient to write the LLG equation in a di- 
mensionless form by defining m = M/M Sl scaled field 
tit = Ht/M s = h e ff + h, scaled time t = t7(|7|M" s ), 
and w(m,h) = W (M , H) / '(^qM 2 ) . The equation for m 
becomes 

(1 + a2 )~~^~ — ~™ x ti — oafi x (m x h t ). (3) 

As shown in Fig.^,, the first term in the right hand side 
of Eq. describes the precession motion around the to- 
tal field and the second term decides the dissipation di- 
rection toward the total field. In our analysis, we always 
use parameters of Co as our references. For Co film|ll|. 
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M s = 1.36 x 10 6 A/m, thus the time unit is approximately 
3.33ps and the energy unit is 2.32 x 10 6 J/m 3 . 

In terms of 8 — 0, the dynamics of the magnetiza- 
tion is determined by the following non-conservative 
two-dimensional nonlinear autonomous dynamical 
equations 0, 



by[K 



(1 + a 2 ) 
(1 + a 2 ) 



-a 



dw 1 dw 
89 ~ iin~fl d^' 
1 dw a dw 



sin 9 89 sin 2 9 



(4) 



Different particles is characterized by different magnetic 
anisotropic energy density functions w{fh,h). In our 
analysis, we always assume it to be the uniaxial with 
its easy axis along the x-direction. Due to the rotation 
symmetry around the easy axis, the applied field can be 
chosen in the x-z plane without losing generality. The 
general form of w(fn, h) can be written as 



w(m, h) = f(m x ) - m x h x - m z h z 



(5) 



where h x and h z are the applied magnetic field along x- 
and z-axis, respectively. 

This naive-looking nonlinear dynamical equation does 
not have exact solutions yet. Analytical results can only 
be obtained in some special cases. For example, analyt- 
ical solutions were foundjU Hjj in the absence of the 
internal effective field (/i e // = 0). The analytical anal- 
ysis can also be done if there is no energy dissipation 
as what is shown in Reference 0,0]. In general, one 
can use the fourth-order Runge-Kutta method to study 
the system numerically. In this study, we will use follow- 
ing f{m x ) in Eq. JSJ with different k 2 and k 4l whenever 
we need to use numerical results to illustrate our under- 
standings, 



f(m x ) = -ljk 2 m 2 x - ^k 4 m\. 



(6) 



Here k 2 , k 4 , accounting for the strength of the anisotropy, 
are positive numbers. In the next section, we shall 
present our main findings. 



III. RESULTS AND DISCUSSIONS 



dw 
dm x 
d 2 w 
dm 2 . 



-0, 
= 0, 



(7) 
(8) 



with 



m~ + m z — 1. For a widely studied case of k 2 7^ 
and k 4 = 0, the S W- limit [Hj is 



(WM 2/3 + (h z /k 2 f/ 3 



(9) 



corresponding to the solid-line in Fig. |3 The SW-limits 
for various choices of k 4 are also plotted in Fig. [21 
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FIG. 2: The SW-limit for various choices of fe, k 4 . Solid 
curve: k 4 — 0; dash-dot curve: k 4 = £2/8; dotted curve: 
fe/4; dashed curve: k 4 = ka/2. 
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The original SW-limit was derived in the static 
case As shown in the dynamical Eq. ((3J), the first 
term on the right hand side will lift the magnetization 
away from the xz-plane. Thus the assumption of the 
SW-limit that the magnetization moves in the xz-plane is 
only true when this term can be neglected. This will hap- 
pen when damping constant becomes infinite (a — > 00). 
In this case, the magnetization will move toward the to- 
tal field as denoted by the big circle passing through the 
north-south poles in Fig. This is the steepest en- 
ergy descent path for the magnetization. Thus, the min- 
imal switching field of SW-limit corresponds to the one 
at which there is only one minimum in energy landscape. 



III.l The exactness of SW-limit at infinitely large 
dissipation 



III. 2 Critical value of damping constant 



The conventional method is based on damping mech- 
anism. Its classical result is the so-called SW-limit. For 
an uniaxial model with the easy axis along x-axis and 
magnetic field in xz-plane, the SW-limit is obtained by 
assuming that the magnetization moves in the xz-plane 
during its reversal. The minimal switching field is given 



In a realistic system, as damping constant is not in- 
finitely large, the magnetization does not need to move 
along the steepest energy descent path. As a result, a 
system may still move from the initial state to the lo- 
cal minimum located near the target state even when 
an external field is smaller than the SW-limit. So, after 
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the external field is removed, the system will move to- 
ward the target state through a ringing mode, achieving 
the magnetization switching. As it was shown in many 
previous studies[H[3, the minimal switching field can 
be smaller than the SW-limit. Numerical calculations [l7| 
show that when the damping constant a < 1, magnetiza- 
tion switching can occur well below the SW-limit. While 
a > 1, the minimal switching field is the SW-limit. Thus, 
it implies a critical a c exists, above which the minimal 
switching field is given by the SW-limit. In Reference 17, 
a c = 1. To the best of our knowledge, there is no clear 
understanding of this result. In fact, previous studies 
seem suggesting that a c — 1 is very special. We want 
to show that there is indeed a critical damping constant 
for a given magnetic anisotropy. But, this critical value 
can be different for different anisotropy, and a — 1 is not 
special. We shall provide an explanation to the a c . 

In order to understand the origin of the critical a, let us 
consider energy landscape under different external field. 
As we mentioned in the previous section, there is only 
one stable fixed point when h > hsw- Asymptotically, 
the system shall always end up at the fixed point for any 
non-zero damping. Thus, if one switch off the field after 
it reaches the fixed point, the system will surely move to 
the target state (state B). There is also a hi < hsw at 
which the energy of system at the initial state equals that 
at the saddle point between two stable fixed points. Thus 
there is no way to switch the magnetization when h <h\ 
because the energy of the initial state (A) is too smaller 
to overcome the potential barrier between the initial and 
final states, hi can be determined from the following 
equations 

0, (10) 



dm x 
w(m x ) = wa, 



(11) 



where wa is the energy at the initial state (m x = 1). For 
a field h between hi and hsw, hi < h < hsw, there exist 
two stable fixed points with a saddle point in between. 
Furthermore, the energy of the initial state is higher than 
that of the saddle point. Fig. [21 is a schematic 3D plot 
of the energy landscape for the case of hi < h < hsw- 
Point A denotes the initial state, whose energy is sup- 
posed to be higher than that of the saddle point (SP) 
between two local minima. In this case, the flow start- 
ing from A will finally evolute into either of two fixed 
points, depending on the value of a. When a is infin- 
ity, the system will evolute into the minimum near the 
initial state along the steepest descent path, as shown 
by line Rl. For the opposite extreme of zero damping 
(a = 0), the system will move along an equal poten- 
tial contour (line R4) surround the two minima (fixed 
points), as investigated in Reference 19. For small a, the 
magnetization can make many turns around the two lo- 
cal minima before it falls into either one. So there is a 
special a = a,; with which the system just touches the 



saddle point (SP) when it rolls down from A, denoted by 
dotted line R3. For a > a,, energy damping is too strong 
for the system to 'climb' over the saddle point, denoted 
by line R2. Value depends obviously on the magnetic 
field, and critical damping constant a c is the value of a, 
at h = hsw- 




FIG. 3: (Color online) The schematic 3D energy landscape 
plot of the case hi < h < hsw- Point A denotes the initial 
state, whose energy is supposed to be higher than that of the 
saddle point (SP). Lines Rl, R2, R3 and R4 show schemati- 
cally four typical evolution trajectories for a — oo, > an, an, 
and 0, respectively. 



One may also understand the result from Fig. 0] of 
trajectories of various a in the energy contour plot at 
h = hsw along 135° to +x-axis. The result in the fig- 
ure is for the uniaxial model with k 2 = 2 and k 4 = 0. 
As mentioned early, the saddle point and one minimum 
merge together at hsw to form an inflexion point de- 
noted by T in Fig.^J It is clear that all trajectories with 
a > a c passing through the saddle-inflexion point while 
all those with a < a c do not. One may notice that all 
curves of a > a c do not move after they reach point T 
because T is a saddle point. But any small fluctuation 
will result in the system to leave T and to end up in FP. 

In order to demonstrate the correctness of our rea- 
soning of existence of a c , and the value of the critical 
damping constant varies with the magnetic anisotropy, 
we carried out numerical calculations on the uniaxial 
magnetic anisotropy model with different ratio of k^jki- 
Fig. El is a-dependence of the minimal switching field for 
ki/k 2 — 0; 1/8; 1/4; 1/2, respectively. Indeed, all curves 
(depend only on the ratio of fci/fe) saturate to their cor- 
responding SW-limit values hsw after a is greater than 
certain values, the critical damping constants a c . Fur- 
thermore, a c is different for different fci, varying from 
a c = 1 for k 4 = to a c = 0.94 for fc 4 = fc 2 /2. Thus, 
a c = 1 is not special at all! 
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FIG. 4: (Color online) The contour plot of w((f>,8) at h = 
hsw for the uniaxial magnetic anisotropy model with = 2 
and &4 = 0. The field is along 37r/4 to +x-axis. Point A is the 
initial point. FP denotes the stable fixed point. T denotes 
the inflexion point. All flow trajectories of a > a c touch T 
while those with a < a c do not. 
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FIG. 5: The minimal switching field versus damping constant 
a. The field is along 3n/4 to +x axis. Square: k^jk^ = 0; 
Circle: k^jki — 1/8; Down-triangle: A^/fe = 1/4; Upper- 
triangle: ki/k2 = 1/2. Smooth connection curves are to guide 
eyes. 

III. 3 Changes of basins under an external field 

We would like to investigate numerically how fixed 
points and their basins change under a field along dif- 
ferent directions. In the following calculation, k% = 2, 
k± = 0, and a = 0.1 is used. In the absence of an ex- 
ternal magnetic field, there are two fixed points at (9 = 
7r/2, 4> = 0) and (9 = ir/2,4> — fl"), respectively. They 
are also the initial state (A) and the target state (B). 
Depending on the direction and strength of an applied 
field, two fixed points, FP1 and FP2, may move away 
from A and B. The numerical procedures are as follows: 
Divide the 9 - (f> plane (9 e (0,tt), and 4> £ (-tt/2, 3tt/2)) 
into 100 x 50 meshes. Each mesh site represents one par- 
ticular state. Starting from an arbitrary state, the state 



belongs to basin of FP1 (FP2) if the system evolutes into 
fixed point FP1 (FP2) after a long enough time. Basin 
FP1 (FP2) is colored white (grey). Fig. HjFig. |S]are the 
basins with different external fields. 
A. Field along the easy axis 

Fig. [fj] (a) is the basin without an external field. The 
basins divide the phase plane into two equal-area parts 
with the separatrix lines of <fi = tt/2 and 37r/2. When 
a field parallel to the easy axis is applied, Points A and 
B are always two fixed points of the system. However 
the areas of their basins vary. Fig. (b)-(d) show how 
the area of basin A shrinks while that of B expands as 
the field increases. When the field reaches the SW-limit 
(Fig. Hjji), the area of basin A shrinks to zero. Thus B 
becomes the only stable fixed point and its basin is the 
whole phase plane when the field is above the SW-limit. 
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FIG. 6: The basins for &2 = 2, ki = 0, and a = 0.1 with 
various fields parallel to the easy axis. Points A and B, which 
are also the fixed points, are the initial and the target states. 
The white color is the basin of fixed point A and the grey 
color is that of B. (a) In the absence of an external field; (b) 
h x = -0.5 ; (c) h x = -1.5; (e) h x = -2 (SW-limit). 

B. Field perpendicular to the easy axis 

When a field perpendicular to the easy axis, say z- 
direction, is applied, the variations of the fixed points 
and basins are shown in Fig. [TJ The positions of two fixed 
points shift along the #-axis, symmetrically located on 
the two sides of +z-axis. It is evident from Fig.[3(a)-(c) 
that the basin shapes become layer-like structures. The 
outer parts of the basins becomes layers first. As the field 
increases, more basin areas become layers. However, the 
total areas of two basins are equal due to the symmetry 
of a perpendicular field when the field is below the SW- 
limit. When the field is higher than a certain value, point 
A will fall into the layered region (Fig. 0:) and may flow 
toward the right hand side of the phase plane. Once the 
magnetization passes the separatrix line (0 = 7r/2), it 
will evolute into the target state B if one can switch off 
the field. This precessional procedure was employed in 
experiments |l l) and numerical calculations |13(. 
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It is interesting to notice that the width of these layers 
become generally thinner when they are away from their 
cores (central parts of basins around fixed points), and 
the number of layers increases with the field. The width 
of a layer is determined by the size of basin core, energy 
variation in the phase plane, and the energy dissipation 
during a 360°-precession (<j> moves from — 7r/2 to 37r/2). 
As the field increases, the areas of the core basins shrink, 
and the precession period is short as well. Thus layers 
become thinner while the layer number increases. The 
energy dissipating rate is ^ = — 1H " a ^ |m x h t \ 2 . The 

precession period can be estimated as about (m x h t )~ l 
since rfi x h t is the angular frequency along the latitudes. 
So the energy loss is proportional to \rh X ht\ — \h t sin#|, 
which is small for 9 is near ir. Thus the width of the 
outer layers is thinner than that of inner ones. Fig. 01 
shows when the field strength reaches SW-limit, the two 
fixed points merge to the same point 9 = (north-pole). 



limit (Fig. 031), the area of the basin of FP1 shrinks to 
zero. 




FIG. 7: The basins in the 9— <j> plane with different strength of 
a field perpendicular to the easy axis, say z-direction. Points 
A and B denote the initial and the target states. The white 
color is the basin for fixed point FP1 and the grey color is 
that for FP2. (a) h z = 0.1; (b) h z = 1; (c) h z = 1.5; (d) 
h z = 2 (SW-limit). Other parameters are the same as Fig.HJ 



C. Field at 135° to the easy axis 

We now investigate how the basins change under a field 
with 135° to the +x-dircction since the minimal switch- 
ing field is smallest around this direction |l7j. Fig. is 
the plot of fixed points and their basins at different field 
strengths. It has both features of cases with parallel 
and perpendicular fields. Increasing the field strength, 
the fixed points shall shift along 0-axis while cf> does not 
change. Unlike the case of perpendicular field, the two 
fixed points have different values. One can see that 
basin FP1 shrinks while that of FP2 expands as the field 
strength increases, a feature with parallel field. How- 
ever, the layer structure also occurs in the outer parts of 
the basins. The stripes become thinner when the field 
increases and 9 is away from the fixed point. At the SW- 




FIG. 8: The basins in 9 — (j> plane with different h along 135° 
to the +X axis. Points A and B denote the initial and the 
target states. The white color is the basin for fixed point FP1 
and the grey color is that for FP2. (a) h = 0.1; (b) h = 0.5; 
(c) h = 0.9; (d) h = 1 (SW-limit). Other parameters are 
the same as Fig. |5] It demonstrates that A and B could be 
connected by one phase flow trajectory in a field with both 
non-zero x- and z-components. 



III.4 Ballistic path 

After seeing the change of fixed points and their basins 
under different field in the previous section, let us in- 
vestigate the field and the reversal time for the ballistic 
connection. Without dissipation, LLG equation is a con- 
servative system. A phase flow is an equipotential curve. 
As it was pointed out in Reference [19j, only a perpen- 
dicular field is possible to connect the initial and the 
target states ballistically. Different from the conserva- 
tive case^lj, the system starting from A will never pass 
through the target state B in the presence of dissipation. 
Even under an infinitely large field, the energy loss dur- 
ing a 180° precession is not negligible. Although 180° 
precession time r decreases as inverse of magnetic field, 
r ~ 7r(l + a 2 )/h when h >> 1 and a << 1, the energy 
dissipation rate goes as ^ = — 1 _ | f a2 \rh x h t \ 2 cx h 2 , thus 



the energy loss during r is proportional to h field|24j|! In 
order to connect both A and B ballistically, one has to 
create a small energy difference between A and B such 
that the energy dissipated on its way from A to B equals 
the energy difference. 

On the other hand, Eq. J2J can be solved exactly 
in the absence of magnetic anisotropy (&2 = 0, = 
0)[13| with solution <j> = ht/(l + a 2 ) and cos6> = [(1 + 
cos o ) e ^ht/(i+ a 2 ) _ ! + cos o ]/[(i + cos 9 )e 2aht ^ 1+a ^ + 
1 — cos 9q] , where 9o is the initial angle between the field 
and the magnetization (we assume the field is along the 
z-axis and the initial <f> is zero). Thus the ballistic reversal 
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corresponding to apply a field along direction 9 satisfying 
-cos 6» = [(1 + cos6»)e 2Q7r - f + cos0]/[(l + cos6>)e 2Q7r + 
1 — cos 9} . It is interesting to notice that the solution is 
unique and the angle is tan(0/2) = e Q7r / 2 . 

Given a damping constant a and magnetic anisotropy, 
described by k 2 and k±, the 180° precession time r(h,(3) 
is a function of field strength h = y/h 2 . + h 2 and its angle 
[3 to the z-axis (J3 relates to 9 by 9 = ir/2 + (3). Thus the 
energy dissipated Ae(h, (3) — Jq ^§dt during r is also a 
function of h and f3. For h » 1, one can neglect the 
magnetic anisotropy, and above isotropic solution should 
be good. Under the limit, Ae is 2/i[l/(l + tan 2 (/3/2 + 
tt/4) e~ 2Q7r ) - 1 / ( 1 + tan 2 (/3/2 + tt/4) )] . The energy differ- 
ence AE(h,(3) between A and B is 2/isin( / 9). Therefore, 
a ballistic path must satisfy Ae = AE (a necessary con- 
dition but not a sufficient one). Due to the symmetric 
reason, one needs to consider only (3 £ (0,7r/2). With- 
out energy dissipation, the only solution is j3 = and 
any h larger than certain minimal value. With large field 
(h >> I) and energy dissipation (a ^ 0), the approx- 
imate solution is tan(/3/2 + tt/4) = e a71 / 2 , the same as 
the isotropic solution tan(6>/2) = e Q7r / 2 . For a / and 
^2 ^ 0, we cannot solve Ae = AE analytically. The field 
configuration of the ballistic connection between A and 
B was found numerically. The results were displayed in 
Figure El 

Surprisingly, the field can be applied in a range of di- 
rection, i.e. a direction window. Given a [3 in this direc- 
tion window, h is uniquely determined. Both the lower 
and the upper bound of this /3-window increase with the 
damping constant. Figure^ is the plot of the upper and 
the lower bounds of f3 as a function of a. The solid line 
is tan(/3/2 + 7r/4) = e air / 2 . Indeed, the window is around 
the solid line. The width of the window depends both on 
the damping constant and the magnetic anisotropy. At 
the zero and the infinite damping constant, the width is 
zero. It is also zero in the absence of magnetic anisotropy 
as illustrate by the exact solution mentioned in the early 
paragraph. Thus, the width is expected to oscillate with 
a for a given magnetic anisotropy. This oscillation was 
indeed observed in numerical calculations as shown in 
Fig. Eb for k 2 — 2, fc 4 = 0. The upper-left inset of Fig. 

is the field and the corresponding reversal time in the 
direction window for a — 0.1 and k 2 = 2, k^ — 0. In this 
particular case, (3 is between 0.134 and 0.156. One sees 
that h increases while reversal time decrease with (3. The 
similar plot for a = 1 is shown in the lower-right inset 
of Fig. EJd. Opposite to the case of small a(= 0.1), h 
decreases and reversal time increases with (3. Thus one 
should compare the lower bound for a < 0.57 and the up- 
per bound for a > 0.57 with tan(/3/2+7r/4) = e a */ 2 since 
it is expected to be exact for h — > oo when the magnetic 
anisotropy can be neglected. An excellent agreement was 
shown in Fig. Et- Fig. is the window width A/3 
as a function of a. Our numerical results indicate that 
the perpendicular configuration employed in the current 



experiments 0,^2 cannot achieve a fully ballistic rever- 
sal. It should be pointed out that above results are for 
the precise ballistic magnetization reversal. As we men- 
tioned early, other field can also switch magnetization if 
one will also like to use the ringing effect at certain stages 
during the reversal process. 



P 



1.5- 








(a) 




1.0- 








y^tan(p/2+K/4)=e a,l/2 


0.5- 






0.0- 







0.0 
0.10- 



0.5 1.0 

a 



1.5 



0.08 

CO. 

< 

£ 0.06 
g 0.04 

o 

"O 

c 

■> 0.02 



0.00 




(b) 



o\ , — , — , — , — lo 

0.13 0.14 „ 0.15 0.16 





a=l 








t 




P 





1 16 1 1R 1 ?0 1?,? 1 74 I 



0.0 



0.5 1.0 

a 



1.5 



FIG. 9: (a) The upper and lower bounds of j3 as a function of 
damping constant a. The solid line is tan(/3/2 + 7r/4) = e Q7r ^ 2 , 
which should be compared with the upper bound of /3 for 
a < 0.57 and the lower bound for a > 0.57. The magnetic 
anisotropy is = 2 and &4 = 0. (b) The window width A/3 
vs. a for = 2 and — 0. Insets: the magnetic field and the 
corresponding reversal time as a function of (3 in the ballistic 
direction window. [3 £ (0.134,0.156) for a — 0.1 (upper-left 
inset), (3 G (1.165, 1.243) for a = 1 (lower-right inset). The 
dashed lines are used to guide eyes. 



III. 5 Discussions 

The SW theory is a classical workjlOj which is based 
on the energy consideration. Its connection to the LLG 
theory was not discussed in literatures. We clarify that 
the SW theory is the infinite dissipation limit of the 
LLG. Under the limit, the magnetization moves along 
the steepest energy descent path. The discrepancy be- 
tween the SW-limit and the numerical minimal switching 
field has been known for many years (l7^. Also, the exis- 
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tence of a critical value of a was known numerically 1 1 7j . 
But its ori gina l was n ot g iven in the previous theoretical 
studies [Hill IT7I HI Hif . To the best of our knowledge, 
its meaning is first revealed now. 

So far we have reformulated the magnetization reversal 
in terms of the language of nonlinear dynamics. The con- 
ventional reversal technique is to make the target state to 
be the only fixed point while the fast switching method 
is to connect the initial and the target states on the 
same phase flow trajectory. Dissipation and magnetic 
anisotropy play very interesting roles in a ballistic mag- 
netization reversal. Without dissipation, only perpendic- 
ular fields larger than a minimal switching field can create 
a ballistic path between the initial and the target states. 
In the absence of a magnetic anisotropy and «/0, field 
direction is not perpendicular but unique, and the field 
magnitude is arbitrary. In the presence of both dissipa- 
tion and anisotropy, allowed directions for the ballistic 
connection form a window. But inside the window, the 
field magnitude is single valued. Due to the energy dis- 
sipation, applying a perpendicular field to the easy axis 
of a uniaxial Stoner particle cannot directly connect ini- 
tial and final states by a phase flow trajectory. Thus, 
the configuration employed in the current fast magneti- 
zation reversal experiments^], ^2 * s n °t the best one. A 
proper field along the direction window can improve the 
best experimental numbers up to now. 

There are advantages and disadvantages for the con- 
ventional and the fast reversal methods. Although the 
switching time in the conventional technique is at the 
nanoseconds, the technology is less demanding and easy 
to implement. On the other hand, the fast switching 
method could make the switching time at picoseconds 
scale, but it needs to have precise control of the mag- 
netic pulse. If it is implemented, the cost should be much 
higher than that of the conventional one. Besides time 
issue, minimal switching field is another concern. The 
minimal switching field in the conventional reversal tech- 
nique occurs when only one energy minimum exists in 
the basin of the final state. Depending on the dissipation 
rate, the minimal field in the precessional fast switch- 
ing method may not be much smaller than the SW-limit. 
So far, in both conventional and fast magnetization re- 
versal schemes, both the direction and the magnitude of 
the magnetic field during the pulse duration are assumed 
to be fixed. If the direction of the magnetic field is al- 
lowed to vary with time, the minimal switching field can 
be even smaller. However it will be more demanding to 
implement technologically. 

IV. CONCLUSION 

In conclusion, we show that the magnetization reversal 
can be conveniently examined in the terminology of non- 
linear dynamics. The presence of the dissipation will not 



fail the fast switching method. We clarify that the SW 
result is the limit of LLG theory with an infinitely large 
dissipation. In this limit, the magnetization moves along 
the steepest energy descent path. We show that there 
is a critical value of the damping constant for a given 
magnetic anisotropy, above which, the minimal switching 
field is the same as the prediction of SW theory. Based 
on the change of fixed points and their basins under an 
external field in different directions, we show that the 
magnetization reversal time can be much shorter when 
the field drives both initial and target states away from 
fixed points, but puts them on the same phase flow trajec- 
tory. In the absence of energy dissipation, the field should 
be applied perpendicularly to the easy axis in order to 
achieve a ballistic magnetization reversal. However, with 
both energy dissipation and a magnetic anisotropy, the 
field can be applied along a direction window. The width 
of the window depends on both the damping constant and 
the magnetic anisotropy. It is zero for either zero damp- 
ing constant or zero magnetic anisotropy (&2 = and 
&4 = 0). Unlike the conventional magnetization reversal 
method, the new scheme demands a precise control of 
picoseconds pulse of a magnetic field. 
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